The third-order law for increments in magnetohydrodynamic turbulence with 

constant shear 
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We extend the theory for third-order structure functions in homogeneous incompressible magne- 
tohydrodynamic (MHD) turbulence to the case in which a constant velocity shear is present. A 
generalization is found of the usual relation [Politano and Pouquet, Phys. Rev. E, 57 21 (1998)] 
between third-order structure functions and the dissipation rate in steady inertial range turbulence, 
in which the shear plays a crucial role. In particular, the presence of shear leads to a third-order 
law which is not simply proportional to the relative separation. Possible implications for laboratory 
and space plasmas are discussed. 
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I. INTRODUCTION 

A well known result in hydrodynamic turbulence the- 
ory is the Kolmogorov-Yaglom ("4/5") law that relates 
the third-order structure function to the energy dissipa- 
tion rate [J 0, 01 ■ Often regarded as a rigorous result 
of the fluid equations, this law requires assumptions of 
isotropy, homogeneity, and time stationarity of the statis- 
tics of velocity increments 5u = u(x + r) — u(x) (veloc- 
ity u, spatial positions x + r and x). In addition, and 
crucially, it also requires adoption of the von Karman 
hypothesis [i| that the rate of energy dissipation e ap- 
proaches a constant nonzero value as Reynolds number 
tends to infinity. Without the need for assuming isotropy, 
one finds 



d_ 

dr, 



{S Ul \6u\ 2 } = -4e, 



(1) 



where (• • ■ ) indicates an ensemble average and a sum on 
repeated indices is implied. If isotropy is further assumed 
then, 



(<5ul|<5u| : 



(2) 



where d is the number of spatial dimensions and Sul = 
f • Su is the increment component measured in the direc- 
tion of the unit vector f parallel to the relative separation 
r. Extension of the third-order law to the case of incom- 
pressible MHD was reported by Politano and Pouquet 
[5j, who remained close to the approximations made in 
the hydrodynamic case. Without assuming isotropy, they 
found 



9 IK 



±|2 



) = -4e d 



which, after adoption of isotropy, reduces to, 
(SzZ\5z±\ 2 ) = -~e ± r, 



(3) 



(4) 



where Sz^ = z ± (x + r) — z^ (x) are the increments of the 
Elsasscr variables and Szf = f • <5z ± . The constants 



are the mean energy dissipation rates of the correspond- 
ing variables z ± = u ± b, where b is the magnetic field 
fluctuation in Alfvcn speed units. 

Here we extend the third-order law in MHD turbu- 
lence to cases in which the isotropy assumption is re- 
laxed. This is accomplished by introducing homogeneous 
shear in the velocity field, a simplified and well-studied 
approach in hydrodynamics [f| 0, H, Q. In particular, 
it supports departures from strict isotropy and introduc- 
tion of shear without consideration of rigid boundaries. 
MHD third-order laws have been applied to systems that 
may also admit departures from strict uniformity, due to 
coherent large-scale gradients; e.g., plasma confinement 
devices Ul 
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and the solar wind [1J, [13|, |14], |15j . For 
systems like these, the homogeneous shear approximation 
may be a reasonable step towards including such large- 
scale effects in the relevant MHD turbulence scaling laws. 
To this end, our derivation of the MHD third-order law 
will include the effect of homogeneous shear, leading to 
a necessarily anisotropic form for the law. 

More specifically, we find that a uniform shear intro- 
duces new terms in the third-order law, so that one can 
no longer conclude that a particular third-order structure 
function, or even a particular integral of a third-order 
structure function, is proportional to the dissipation rate 
times the relative separation length r. This is in marked 
contrast to the situation for the fully isotropic hydrody- 
namic and MHD cases, given here as Eqs. $2$ and (TJ|. It 
is, however, entirely consistent with the work of Lindborg 
[lq and Casciola et al. [l7[ , who derived modifications to 
the form of the third-order law for hydrodynamics with 
shear. 

The principle theoretical result given below is that a 
uniform shear indeed is responsible for changing the form 
of the third-order law, whereas a mean magnetic field 
does not produce such structural changes. Implications 
for solar wind, laboratory, and astrophysical measure- 
ments of turbulence are suggested, and in particular the 
primacy of the third-order law in unambiguously defining 
an inertial range is challenged. 
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II. ENERGY DECAY WITH LARGE-SCALE 
FIELDS 

The third-order law is often derived from the steady- 
state version of an equation related to energy decay. 
To obtain the version of the law appropriate for MHD 
with uniform velocity shear, we follow the same proce- 
dure used previously for MHD [f| [l3| , combined with the 
method of Casciola et al. [ljj for extending Eq. @ to 
include shear. A uniform magnetic field is also retained, 
although only the simplest of its consequences will enter 
the discussion. 

First, let us employ a Reynolds decomposition of the 
velocity field v = U + u into a mean velocity U(x) and 
a fluctuating component u(x, t), where (v) = U and 
(u) = 0. Similarly we write the total magnetic field, con- 
veniently expressed in Alfven speed units, as B = b+Bo- 
We assume Bo is constant and uniform, but that U(x) 
varies in space. However this variation will be taken as 
non-random and slowly- varying, so that the turbulence 
properties can be treated as locally homogeneous. 

Now we write the incompressible MHD equations at 
two positions, x and x' = x + r: 

d t zf = -{zl + U k T B ok )d k (U t + zf) - d % P + vd k d k zf, 

(5) 

d t zf = -{zf+U' k TB ok )d' k {Ui+zt')-d[P'+vd' k d' k zt'. 

(6) 

Here the prime denotes quantities at position x', P is the 
pressure, and v is the kinematic viscosity, taken equal to 
the resistivity hereafter. Subtracting Eq. (JSJ from Eq. (J6]) 
yields the following equation for the Elsasser increments 
Sz ± =z ± (x') -z ± (x): 

d t 8zf = -(5U k + Sz^d'Jzf 

-(4 + U k T B ok )(d' k + d k )5zf 

~{5zT + SU k )8 k U t 

-{zf + U' k TB ok )S(d k Ui) 

-(flj + di)5P + v{d' k d' k + d k d k )5zf, (7) 

where we use the property that the primed and unprimcd 
coordinates are independent, so that d k zf' = and 
d' k zf = 0. 

As noted above, we seek an equation related to energy 
decay. Multiplying the previous equation by 2Szf 1 and 
averaging yields 

MSzff) = -±.((6U k + SzZ)\6z?\*) 

+{\Szt\ 2 (d k U k + d' k U / k )) 
-2(d k Ui6z?(6z* + SU k )) 
-2{{zf + U l k TB 0k )8{d k U i )8zf) 

+2^(|<Szf| 2 )-M|a fc z±| 2 ). (8) 

In arriving at this expression we make use of d k (») = 
— and d' k (*} = These latter relations fol- 



low from spatial homogeneity (i.e., translation invariancc 
of the statistical properties), which can be considered for 
some systems to be an exact property (see following sec- 
tion) or an approximation, e.g., in the case of a weakly 
inhomogeneous system. The main results here will be for 
strict homogeneity. 

The last term of Eq. ([8|) can be identified with the 
dissipation rates 

e ± = v(\d k z±\ 2 }, (9) 

which for steady state are also the mean energy trans- 
fer rates. Following the usual arguments [J], in the 
limit of vanishing viscosity v — ► 0, it is assumed — not 
proven — that the remain nonzero, and in effect are 
externally prescribed by the rate of supply of turbu- 
lence energy (and cross helicity). Although this non- 
trivial assertion is physically plausible [l8T |. it nonethe- 
less prevents the subsequent developments, including the 
classical 4/5-law, from being considered an exact conse- 
quence of the fluid equations themselves. Furthermore, 
the penultimate term in Eq. ([5]), also involving the vis- 
cosity, is assumed to vanish at high Reynolds number 
when we are examining the inertial range of separations. 
For the above-stated set of approximations, the incre- 
ments r are restricted to lie in the inertial range, that is 
separations smaller than the correlation length (energy- 
containing scale) and bigger than the dissipation scale 
(scale at which fluctuations are critically damped). For 
variations of the set of assumptions that lead to a third- 
order law, see e.g., Hill (l9j |. 



III. MHD WITH HOMOGENEOUS SHEAR 

The above relations need not be strictly homogeneous, 
as variations in U over the slowly-varying large scales 
may be present. To rectify this and arrive at a general 
law that is translation invariant, we now specialize to the 
case of a homogenous shear flow, alluded to in the intro- 
ductory section. With this choice the tensor dUi/dxj is a 
constant matrix independent of position. The turbulence 
is then homogeneous and all terms in Eq. {5]) — both co- 
efficients and averaged terms — are only a function of the 
separation vector r. 

Under the hypothesis of steady-state turbulence, the 
left-hand side of Eq. (|SJ) vanishes. Integrating in r, over 
a volume V that is enclosed by a surface 5, the equation 
becomes: 




= -AVe ± , (10) 

where V is the volume of the region V and n k is a unit 
vector normal to S. 

If the region of integration is a three dimensional 
sphere of radius r, volume V r and surface S r , the in- 
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tcgration yields 



.dUi 



[r k ((8z% + 5U k )\5zf\ 2 )] dS r + 2^ J (Szf5z^)dV r 
~^e\ (11) 



where, in the first term of the equation, fk is the unit vec- 
tor normal to the surface of the sphere, and now in spher- 
ical (r, 8, (f) coordinates dS r = r 2 d(cos6>)d0 = r 2 dfl. 
Equation (|11|) may be interpreted as the integral form 
of the third-order law for incompressible homogeneous 
MHD turbulence with an external velocity field that is 
constant in time but which can vary linearly in space. 
By setting U = and assuming isotropic turbulence, 
Eq. (TTTj) will recover the standard third-order law for 
isotropic MHD turbulence Q, given here as Eq. Ql 

In standard derivations for isotropic turbulence 0, Il3l . 
fl4j |. shear is necessarily lacking, and it is assumed that 
the structure functions are rotationally symmetric. In 
that case the above relation is simplified by carrying out 
the integrals explicitly. (For a more general case, see 
below.) Here we allow for anisotropy induced by cither 
the large-scale magnetic field, or by the imposed homo- 
geneous shear. Note that the large-scale magnetic field 
Bo does not appear explicitly in the third-order relation, 
even though it is well documented that such a field in- 
duces spectral anisotropy in MHD turbulence [20j | . 

We now further specialize to the large-scale homoge- 
neous shear flow U = U x (y)x — ayx in a cartesian 
(x, y, z) system, with a = const. The integral form of 
the third-order relation becomes 



r • <(5z =F | ( 5z ± | 2 )J dn 



ar 3 <j)(r ■ x)(f • y)(\5z ± 
2a J (SztSz*) dV r = - 



dn 



167rr d 



(12) 



Denoting an angular average over a shell of radius r as 
(• ■ • )n arL d a volume average over a sphere of radius r as 
(■ ■ -)v wc may rearrange the above equation as 



((Sz^S^n 



-ar(((f •x)(f-y)| ( 5z ± | 2 ))o 
-Zar{(5z±5zf)) v - ^re±, (13) 

where, again, 5z^ = r-Sz ± . This form, based on a spher- 
ical region of radius r, indicates that all three terms on 
the right hand side of the equation have an explicit pro- 
portionality to r; moreover, the first and second of these 
also admit an implicit dependence on r. The quantity 
on the left side of Eq. ([13]) is the MHD analog of the 
usual third-order structure function that appears in the 
Yaglom and Kolmogorov laws [l], Q , and we see that in 
the presence of homogeneous shear it is not simply pro- 
portional to the dissipation e ± . 

At this point we remark on an alternative form that 
the third-order law can assume that may be revealing 



in anisotropic cases. Recall that Eq. (|10[) is valid for 
an arbitrary volume V and its associated bounding sur- 
face S. The advantage of employing a spherical vol- 
ume V is that when the flux is isotropic, the integrand 
in the surface integral will be independent of the direc- 
tion of r, making the integration trivial. Unfortunately, 
this property is lost when the turbulence is anisotropic 
[20l . l2ll ]. However, provided that the (energy-like) vec- 
tor flux F + = ({6z~ + 5U)\6z + \ 2 ) is smoothly varying in 
r, it is in principle possible to find a set of nested sur- 
faces S(V) [labeled by their enclosing volume V and with 
unit normal vectors rig], such that the normal compo- 
nent of the vector flux F + is uniform across <S(V). Then 
<f s dS hg -F + = F+(V)S, where the constant normal flux 
F n is labeled by the volume V bounded by the surface, 
and S is the value of the surface area. The partner quan- 
tity F is defined analogously. Provided these nested sur- 
faces can be found, the homogeneous shear case, Eq. QXT 
can then be reduced to 



n 



s± ■ ( 5z 
2aV ± 



SXJ)\6: 



,± |2\ 



AV ± 



(14) 



where and are the volumes and associated surface 
areas that admit constant normal fluxes F r ^(V r± ). Note 
that in general the constant flux surfaces S + and S~ are 
expected to be different, e.g., due to cross helicity effects. 

When homogeneous shear is absent the result in 
Eq. (fT4|) reduces to the formal anisotropic third-order 
law 



W ± 



(15) 



The latter can have application in the cases in which 
anisotropy is present due to a mean magnetic field Bo 7^ 
0. 



IV. SUMMARY AND DISCUSSION 

We examined the mixed third-order Elsasser structure 
functions for MHD turbulence, incorporating a constant 
sheared velocity (homogeneous shear) field in addition 
to homogeneous fluctuations, under a set of assumptions 
that parallels those used in standard turbulence theory 
to derive the Kolmogorov 4/5-law. In analogy to the 
findings of Casciola et al. [17| and Lindborg [16( for hy- 
drodynamics, wc find that a law can be obtained for sta- 
tionary homogeneous turbulence that relates third-order 
structure functions and dissipation, but which also in- 
volves additional terms. For MHD with a constant im- 
posed shear, there are shear-related terms that appear 
in this modified third-order law, as in the hydrodynamic 
case. On the other hand, a uniform magnetic field does 
not appear explicitly in this relation. 

On the basis of a very simple estimate we expect the 
new terms in the third-order equation to be of signifi- 
cance when the large-scale velocity increments are of the 
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same order or larger than the fluctuation increments at 
the same separation r, that is, when SU ~ Sz. In some 
applications this condition may be realized, and conse- 
quently the classical third-order law is modified by these 
new terms. We suspect that for solar wind turbulence, 
as well as for laboratory devices, the present generalized 
form of the third-order law will be relevant. In particu- 
lar, the modified MHD third-order law no longer admits 
an interpretation purely in terms of energy transfer and 
dissipation, and therefore differs from the isotropic case 
without shear. 

Further extensions of the third-order law can also be 
undertaken. For example, by including a large-scale but 
non-uniform magnetic field. This will induce additional 
terms in the generalization of Eqs. (fl"TjF(|13ll . 

As a final remark, we note that the modifications of 
the third-order law for energy decay that we describe 
here can be anticipated in the structure of scale-separated 
transport equations derived for MHD in a weakly inho- 
mogeneous medium [22], [23| . These two-scale transport 
equations provide a formalism for evolution of second- 
order correlation functions, and include nonlinear decay, 



analogous to our third-order structure functions, along 
with advection and shear terms. On this basis, one could 
have already concluded that the third-order law requires 
modification in the presence of large-scale shear. The 
present study concentrated only on the special case of ho- 
mogeneous shear, and generalizations of the third-order 
law have been found. 

We expect that future studies based on numerical sim- 
ulations may provide explicit verification and examples 
of the relationships we propose here. Taking into account 
effects like shear, observational studies may prove useful 
in a variety of systems with large-scale shear flows, such 
as astrophysical and laboratory plasmas. 
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